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Abstract 

We consider the three-body decay t\ — > bi> T T + and propose the asymmetry of the r polarization 
perpendicular to the decay plane as a CP sensitive observable. We calculate this asymmetry in the 
Minimal Supersymmetric Standard Model with the parameters /j and A t complex. In the parameter 
domain where the decay ti — > bu T r + is important this asymmetry can go up to ±30%. We also 
estimate the event rates necessary to observe this asymmetry at 90% CL. 



1 Introduction 

The experimental search for supersymmetric (SUSY) particles will have high priority at the upgraded 
Tevatron and at LHC. The precision determination of the SUSY parameters will be the main goal of a 
future e + e~ linear collider JjJ. The analysis of scalar top quarks U, i = 1,2, will be particularly interesting, 
because of the large top Yukawa coupling involved in this system. Due to the effects of the top Yukawa 
coupling the lighter stop may be relatively light if not the lightest charged SUSY particle [^) . 

A phenomenological study of production and decays of 3rd generation sfermions at an e + e~ linear 
collider with cms energy in the range -^i = 0.5 — 1 TeV has been given in Ref. ||. There it has been 
shown that by measuring production cross sections with polarized beams the masses of the top squarks 
mj , TOj 2 and their mixing angle 0^ can be determined. It has also been shown that with an integrated 
luminosity of 500/fe -1 a precision of about 1% can be achieved. The precision to be expected for the 
underlying SUSY parameters has also been estimated. This analysis has been performed in the Minimal 
Supersymmetric Standard Model (MSSM) with real parameters. 

However, the assumption that all MSSM parameters are real may be too restrictive. In principle, the 
Higgs-higgsino mass parameter /i and the trilinear scalar coupling parameters Af of the sfermions / may 
be complex. These complex parameters are new sources of CP violation and may provide potentially 
large SUSY contributions to the electric dipole moments (EDM) of electron and neutron. The very small 
experimental upper limits of the electron and neutron EDMs, therefore, may lead to restrictions on the 
complex phases. Recent analyses have shown that in mSUGRA-type models the phase of /i is restricted to 
\<Pu\ % 0.1 — 0.2 for a universal scalar mass parameter M < 400 GeV, whereas the phase of the universal 
trilinear scalar coupling parameter A is correlated with (p^, but otherwise unrestricted Q. One can 
conclude that in models with more general parameter specifications the phases of the parameters Af 
of the 3rd generation sfermions are not restricted at one-loop level by the electron and neutron EDMs. 
There may be restrictions at two-loop level || . Furthermore, a complex trilinear coupling parameter A t 
in the stop system can also lead to interesting CP violating effects in top quark production, as discussed 
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For a complete analysis of the stop sector one has to take into account that the parameters /z and 
A t may be complex. The parameter \fj,\ and its phase ip^ will presumably be determined by measuring 
suitable observables of the chargino and neutralino sector [j7j. For the determination of \At\ and its phase 
ipA t appropriate observables in the stop sector have to be measured. However, it may be difficult to define 
a suitable CP sensitive observable in stop decays if the main decay modes are two-body decays. 

In the present paper we define a CP sensitive asymmetry in the decay t\ — > bv T r + . As this is a 
three-body decay the polarization of the t + normal to the decay plane is sensitive to CP violation. The 
appropriate CP sensitive observable is defined by the asymmetry of the r polarization perpendicular to 
the decay plane. As we will show this asymmetry can go up to 30%. Moreover, we show the existence 
of parameter regions where the decay t\ —> bv T T + has a sufficient branching ratio allowing for the 
measurement of this asymmetry. We perform our analysis in the MSSM with fi and A t complex. We 
focus on scenarios where only the decays t\ — ► bi>t£, t\ — ► bivg |@|, t\ — > c\\ Q and four body decays 
JTc| are kinematically allowed. We assume that the lightest neutralino Xi is the lightest SUSY particle 
(LSP). The CP asymmetry defined above is analogous to that defined in Ref. [jll) in case of top quark 
decays. 

In Section 2 we shortly review stop mixing in the presence of complex parameters. In Section 3 we 
give the formulae of the CP violating observable. In Section 4 we present numerical results for the phase 
dependences of the CP asymmetry. We give an theoretical estimate of the event rates necessary to observe 
the CP sensitive asymmetry at an e + e~ linear collider with s/s = 0.5 - 1 TeV. Section 5 contains a short 
summary. 



2 t L - t R Mixing 

We first give a short account of tL - in. mixing in the case the parameters /z and A t are complex. The 
masses and couplings of the i-squarks follow from the hermitian 2x2 mass matrix which in the basis 
(tL,tfi) reads 

/ M? ll e-^\M?J \fh\ 

Cm = -(*!, t R )\ > (!) 

V e***\M? | Mf \i R 

V 1 iLFt' IBS ' V / 

where 

M? ll = M| + (i-^sin 2 e w )cos2/3m|+m?, (2) 

M? rr = + ^ sin 2 ® w cos 2/3 m% + m^, (3) 

M? rl = (Ml R y = m t (A t -n*cot(3), (4) 

ft = ar g[^t - V* cot 0\ , (5) 

where tan/3 = 1)2/1)1 with 1)1(1)2) being the vacuum expectation value of the Higgs field H^H®), m t is 
the mass of the top quark and <3>w is the weak mixing angle, /i is the Higgs-higgsino mass parameter 
and Mq, Mfy, A t are the soft SUSY-breaking parameters of the stop system. The mass eigenstates t$ are 

(tiM) = (i L J R )n iT with 

e llfi cos 61 sin % 
R'=\ ' J. ((,) 

— sin Of e~ ZVi cos 9 
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with 

-\M? I M? -ml 



(7) 



* J\M? | 2 + (m? - A/? ) 2 ' ' * ,/|M? | 2 + (m? - M? ) 2 ' 

V 'LH tl ILL, V t-L-R f l ILL 

The mass eigenvalues are 



ii, — in mixing is naturally large because of the large top quark mass entering in the off-diagonal elements 
of the mass matrix (see Eqs. (|l|) and (jl])). This is important for the CP sensitive observable discussed 
below, because it is proportional to sin 0j cos Qi . Note further that for \At \ 3> \n\ cot f3 we have iff w 



3 Tau Polarization Asymmetry 

The parts of the Lagrangian relevant for the three-body decay t\ — > bv T r + are 

= <?&(M>£ + + h.c, (9) 

£t*,* t = S*f {kf*PR + lj r *Pz)r K + h.c, (10) 

where 

l\. = -e*Pi cos%^i + F t sin%^- 2 , fcy = J^e - *** cos%l7/ a> (11) 

£ T =-V$i, fcf=Y r C/; 2 . (12) 

^ is the weak coupling constant, Pr,l = 1/2(1 ±75) and the Yukawa couplings are Yj = mi/\/2?TJw sin /?, 
Y> = rrib/^mw cos/3, F T = 77i t /\/2tow cos/3. The unitary 2x2 matrices U and V diagonalize the 
chargino mass matrix p^ |. A convenient parametrization of U and V for complex parameters can be 
found in 

The three body decay t\ — > bv T r + proceeds via exchange of charginos xf, i = 1,2. We consider the 
polarization of the outgoing t + perpendicular to the decay plane, which is sensitive to CP violation. We 
define the unit vector 

Pt x Pb h „, 
ejv = j-z ztt- (13) 

\Pt XPb\ 

The average polarization of the outgoing r + in the direction ejy in the decay t\ — > bv T r + is given by 

pT + = flfo ^ bi/ T r+(e Ar )) - gfc -> bi> T T + (-e N )) 
N B(h -» bv T T+(e N )) + B(h -» Myr+(-SW)) ' 

For the CP conjugated process we get 

- _ B(ti -> 6i> t t~ (ejv)) - -* &£ T i~~ (~e N )) 

' N ~ 5 — = • 

_B(<i — > bv T T-(e N )) + — > W t t~{— ejv)) 

Note that Pjy = —VJ^ since the couplings are complex conjugate to each other. The observables Pjy" 
and V T N are odd under naive time reversal TV, where only the polarization and momentum vectors are 
reversed but initial and final states are not interchanged. We can define a CP sensitive asymmetry of 
the form 

Acp = \{V$-V r N ). (16) 
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In order to obtain a CP asymmetry which for practical reasons is more useful, we have to take into 
account the subsequent decay v T — > Xi u t- The complete decay chain is then t\ — > bv T T + — > 6t + Xi^V- 
Another decay chain leading to the same final state is ?i — > &TjZ/ T — » brx\v Tl where in the second step 
the fj decays into rx?- In the fj rest system tj — + t\\ has an isotropic decay distribution. The decay 
mode t\ — > 6t;j/ t — > brx\v r can be easily incorporated in our consideration. 

However, the decay t\ — > 6W^x? ~~ * ^> t x\ v t also leads to the same final state [Q. The decay ?i — » 
6Wx$ — > brxi^r is more involved, because the W polarization leads to a non- vanishing correlation 
between the p T - pb plane and the r polarization. While this may also lead to CP sensitive effects, 
in the present paper we confine ourselves to the discussion of the CP asymmetry in the decay chains 
ti — » bv T T — > brxi^T and t\ — > 67^1/7- — > brxi^r, assuming that < 111^ + m^o + nib- 

We define a CP asymmetry similar to Eq. (|l6|), but for the final state brxl^T- 

1 ,_,r+ 

2 

with 



■Acp = o(^ -^), (17) 



= /(fly)) - B(h -+ f(-e N )) V , T - = gfc -> /(eW)) - gfc -> /(-eW)) 

W " Bit, ^ f(e N )) + B(h ^ f(-e N )y 

where we have introduced a shorthand notation /(e/v) = Xi^rb r + (eV) and /(ejy) = Xi^rb t~ (ejv)- The 
branching ratio for the decay ti — ► Xi^rb r is to a good approximation given by 

2 

-> X> T & r) ~ B(ti -> bv T r)B{v T -» x> T ) + ]T B(ti -» br^B^ -> (19) 

This can also be seen by using the formulas for the 4-body stop decays given in |L0| . Relation (|l9| ) holds 
if Irrifj — rnv, T \ ^> + r^, which is naturally fulfilled. If B(v T — ► Xi^r) — 1 and B(f$ — > Xv^) — lj 
then Eq.jl^) can be rewritten as 



, r(fi -» 6i> r T + (eAf)) - T(fi -» fa> r T + (-ejv)) - -» &t> T r (e N )) + T(ti — > &£> t t (-e N )) 
■^cp — TT^. ! 

^ r unpol 

(20) 

where r unpo i = Y(t\ — ► 6t> T r) + X^=i r(ti — > bfiV T ). A very useful approximation for A' CP can be 
obtained in the limit (p^ — > 

-Acp - : pr^ 2 m ?i r t r T | sin2%| sin^- T, (21) 

-L unpol 



where 



1 f d 3 p b d 3 p T d 3 pp T 4 |p T ||p b |sin6» feT 

(2tt) Sip^-pb-pr-p^j-, - r -, (22) 



2m fl 7 (271)32^(271)32^(2^2^/ ' ^ ^ " — (jA ± - mi ± )(pi ± - mU) 

where p^ = (jb^ , 0) , p^± — pi x —pb, m ^± are the masses of the charginos and sin 6b T is the angle between 

the b quark and the t lepton. Eq. (^TJ) is very instructive, because it exhibits the main features of the 
behaviour of A' CP . For example, the sign of A 'g P is given by sgn(A' CP ) = sgn(ifi). The unpolarized 
decay width r unpo i appearing in Eq. ( |2(i| ) and ( pl| ) is given in [|| for real parameters and in Q for 
complex parameters. 

Finally, we want to remark that the absorptive part of the amplitude caused by the Breit-Wigner form 
of the chargino propagators is strongly suppressed, because we have assumed that the decay it — > X\b is 
not accessible. Therefore the rate asymmetry is expected to be at most of the order of w 10~ 3 . 
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4 Numerical Results 



In the following we present numerical results for the CP sensitive asymmetry A' CP defined in Eqs. (|T^) and 
(|l8|). Our input parameters are , m ( - 2 , m^ T ,\Af\, (pA t , A T , tan (3, M2, \ assuming Mi = 5/3 tan 2 ®w M2 
with Mi and M2 real. For simplicity we set ip^ = 0. We impose the approximate necessary condition for 
the tree-level vacuum stability \A t \ 2 < i{M^ + M^ + {m 2 A + m 2 z ) cos 2 (3—m 2 z /2) |l6). For the pseudoscalar 
Higgs mass which appears in this condition we choose for definiteness itia = 150 GeV. We have checked 
that in the numerical examples studied below the restrictions from the electron and neutron EDMs at 
two-loop level § are fulfilled. 

First we consider the influence of the parameter ipA t - In Figs.|l]a and [TJb we show A' CP as a function 
of (fA t for Mq > Mjj (| cos%| < I sin#j|) and Mq < Mfj (| cos6y > | sin%|), respectively. We display the 
asymmetry for the four scenarios (\fi\ = 400 GeV, tan (3 = 3), (\p,\ = 400 GeV, tan/3 = 10), = 700 GeV, 
tan/3 = 3) and = 700 GeV, tan/3 = 10), taking m h = 240 GeV, m~ t2 = 800 GeV, mp T = 200 GeV, 
M2 = 350 GeV. We fix the mass of f\ by taking Mg = 0.9M^. As can be seen, the CP asymmetry 
A' CP is much larger for Mq > Mfj than for Mq < Mfj. The reason is, that for Mq < My the t\ has 
a larger component than for Mq > Mfj, implying a stronger coupling to the gaugino component of 
the charginos. This in turn implies a larger r unpo i for Mq < Mfj. Cearly A' CP has to vanish for the 
CP conserving case <*pA t = 0, ±7r (see the corresponding factor in Eq. (pl|) ) . For better understanding of 
the value of A' CP in the region \(fA t \ ^ 7r /4, we show in Fig. || the branching ratio B(ti — > Xi^t^t) as 
a function of ipA t , for the two scenarios (\fj,\ = 400 GeV, tan/3 = 3) and (|^| = 700 GeV, tan/3 = 3), 
keeping the other parameters as in Fig.^ja. The decay width T(ti — » xl^rbr) has a similar behaviour. 
The minimum at ipA t — can be traced back to the fact that l/f-J (see Eq. (|IT|)) has a minimum there, 
due to a negative interference between the gaugino and higgsino contribution. 

In Fig.| we show A' CP for three values of \A t \ = 600 GeV, 1000 GeV and 1300 GeV, assuming 
Mq > Mfj. The other parameters are m~ tl = 240 GeV, m t ~ 2 = 800 GeV, = 200 GeV, A T = 0, 
M 2 = 350 GeV, H = 600 GeV and tan/3 = 3. For \A t \ < 1000 GeV A' CP increases with increasing \A t \, 
because | sin20j| increases (see Eq. (^l|)). The decrease of A' CP for \A t \ > 1000 GeV is explained by the 
fact that r unpo i increases stronger than | sin2#j|. 

The polarization of the r is analysed through its decay distributions. Usually the decay modes 
r — > iris, pv, a\v, iivv, evv are used as analyzers. As we are interested in the transverse polarisation of 
the r lepton, we take only the pv and a\v final states for our analysis. The sum of the branching ratios 
of these two decay modes is « 34%. We take for the sensitivities for measuring the polarization of 
the r lepton the values quoted in Ref. |l7j for an ideal detector. Moreover, the numbers quoted are for 
longitudinal tau polarization and it is expected that the sensitivities for transversely polarized tau leptons 
are somewhat smaller. To account for both effects we assume a reduction of the sensitivity S of 10% [18|. 
Furthermore, we assume that the direction of flight of the r lepton can be reconstructed. Following [17|, 
the error in measuring the polarization asymmetry is given by 

sv ' Tn = T^Tw ' (23) 

where SVod = 0.3, N T — B(ti — > Xi^rbr)^ , with iVj- being the number of produced tit\ pairs in the 
process e + e~ — > txt\. 

The minimum number iVj of produced t\ pairs necessary to observe A' CP at 90% confidence level 
(CL) is then given by 

1 / 1.64 ' 



Nil biI^xIvmAsua'cpV ' (24) 

In Table 1 we display the values of the asymmetry A' CP and the numbers Ni needed to observe this 
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asymmetry at 90% CL. 



Table 1: The CP asymmetry A! CP defined in Eq.([l7|) and the number iVj of t\ti pairs required to 
measure this CP asymmetry at 90% CL, choosing mi — 240 GeV, mi — 800 GeV, m„ T — 200 GeV, 
Mi = 350 GeV and \A t \ = 1000 GeV. For Mq > A% we take = 700 GeV, tan/3 *= 3) and for 
Mq < M G we take = 400 GeV, tan/3 = 10). 



A' 


N h x 10~ 3 


<PA t 






0.05 


31 


tt/2 


Mq 


>M 


0.11 


7 


tt/4 


M Q 


>M 


0.21 


4 


tt/8 


M Q 


>M 


0.033 


48 


tt/2 


M Q 


<M 6 


0.030 


58 


tt/4 


M Q 


<M 


0.018 


165 


tt/8 


Mq 


<M 



The cross section for e + e" — > titi for m^ = 240 GeV is about 5/6 at y/s — 500 GeV and in the 
range 15 - 120/6 at ^/s — 800 GeV, depending on stop mixing and beam polarization. If the branching 
ratio B(ti — > bv T T + ) is of the order of 5% or larger, then there are good prospects to measure the CP 
asymmetry A' CP at a linear collider with an integrated luminosity of 500/6 -1 . 



5 Summary 

In this paper we have studied the CP violating observable defined by the asymmetry of the r polarization 
perpendicular to the decay plane in the three-body decay t\ — > bv T T + . In the parameter domain m^ < 
Tn^± + mb, mi t < mw + m x'l + m b the branching ratio for this decay mode can be quite large. We have 
calculated this CP asymmetry and the branching ratios in the MSSM with the parameters fi and A t 
complex. We give numerical predictions for the case of an e + e~ linear collider with cms energy ^fs = 0.5 
1 TeV. The asymmetry can reach values up to ±30%. We give a theoretical estimate of the number 
of produced titi pairs necessary for measuring this asymmetry at 90% CL. There are good prospects 
to measure this CP asymmetry at a linear collider with an integrated luminosity of 500 fb , if the 
branching ratio of the decay t\ — > bu T T + is about 5% or larger. 
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Figure 1: The CP sensitive asymmetry A' CP as a function of ipA t - The input parameters are 
m h = 240 GeV, m~ t2 = 800 GeV, mp T = 200 GeV, M 2 = 350 GeV, \A t \ = 1000 GeV, (tan/3 = 
3,|/x| = 400 GeV; thick solid line), (tan/3 = 10, |,u| = 400 GeV; thin solid line), (tan/3 = 3, = 
700 GeV; thick dashed line), (tan/3 = 10, = 700 GeV; thin dashed line), for the cases: a)Mg > Mjj 
and b)Mn < Mjy\ 
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Figure 2: The branching ratio of t\ — > Tb\\v T as a function of ifiA t , with Mq > M^m^ — 240 GeV, 
m t ~ 2 = 800 GeV, mz, T = 200 GeV, tan/3 = 3,M 2 = 350 GeV, \A t \ = 1000 GeV for the two cases = 
400 GeV (solid line), |/z| = 700 GeV (dashed line). 




Figure 3: The CP sensitive asymmetry A' CP as a function of fA t - The input parameters are mj = 
240 GeV, m ?2 = 800 GeV,™^ = 200 GeV, M 2 = 350 GeV, tan/3 = 3, |m| = 600 GeV for \A t \ = 
600 GeV (solid line), \A t \ = 1000 GeV (dashed line), \A t \ = 1300 GeV (dotted line). 
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